Delay and distortion of slow light pulses by excitons in ZnO 
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Light pulses propagating through ZnO undergo distortions caused by both bound and free exci- 
tons. Numerous lines of bound excitons dissect the pulse and induce slowing of light around them, 
to the extend dependent on their nature. Exciton-polariton resonances determine the overall pulse 
delay and attenuation. The delay time of the higher-energy edge of a strongly curved light stripe 
approaches 1.6 ns at 3.374 eV with a 0.3 mm propagation length. Modelling the data of cw and 
time-of- flight spectroscopies has enabled us to determine the excitonic parameters, inherent for bulk 
ZnO. We reveal the restrictions on these parameters induced by the light attenuation, as well as a 
discrepancy between the parameters characterizing the surface and internal regions of the crystal. 

PACS numbers: 78.20.e 



The problem of optical pulse propagation in a medium 
has attracted attention from the beginning of the last 
century Loudon [2j studied the transfer of electro- 
magnetic energy in a local dielectric with a resonant ab- 
sorption line. A large number of followers |3|— ISf] extended 
his analysis to include spatial dispersion. The experi- 
mental time-of-flight study of Chu and Wong Q of a 
laser pulse tuned to the bound exciton (BX) line in GaP 
verified the prediction of Garrett and McCumber [7[ on 
the dramatic variation of group velocity at a resonance. 
Their experiments promoted theoretical investi gati ons of 
distortion of the temporal shape of a pulse [8l4lQ|. Re- 
tardation of exciton-polariton propagation was studied 
in various semiconductors [lll-ll3| , possessing an isolated 
resonance, rather than an array of closely situated lines. 

Renovation of the interest in slowed light was stimu- 
lated by efforts in quantum information processing (l^ . 
That requires a delay significantly exceeding the pulse 
duration, while their ratio is limited by ~4 for quantum 
coherent effects [l5| . The resonant optical dispersion al- 
lows one, in principle, to overcome this limit, however 
at the price of severe attenuation of the pulse intensity. 
Wide-gap semiconductors, like ZnO and GaN, promis- 
ing for various optoelectronics applications, demonstrate 
strong excitonic resonances that makes them suitable for 
such a light slowing. Retardation of photoluminescence 
(PL) by 200 ps was observed near the 3.36-eV BX line 
in a 1-mm long ZnO sample [l6|. In GaN, the delay was 
twice higher for the same sample length near BX reso- 
nances ua. In view of these findings, the question about 
the distortion of a slowed pulse has been revived [l8| . 

Here, we report on delay and distortion of light pulses 
by excitons in ZnO. The results are novel for two rea- 
sons: First, ZnO possesses numerous lines related to dif- 
ferent states of donor bound excitons. The fine spec- 



trum of these excitons is probed for the first time by 
time-of-flight spectroscopy via the local distortion of a 
pulse shape at different energies. Second, simulation of 
the general shape of the transmitted pulse permits us 
to determine the parameters of exciton-polariton reso- 
nances, inherent for bulk ZnO. They differ from those 
given by surface-probing techniques. The knowledge of 
correct parameters is important in many aspects; in par- 
ticular, they control the delay and attenuation of a pulse. 

The time-of-flight experiments were performed using 
the pulses of a picosecond tunable laser (Mira-HP, sec- 
ond harmonic). A Hamamatsu streak camera with a 2 
ps temporal resolution was exploited to record the time- 
resolved (TR) images of the pulses. Their temporal 
width (~30 ps) was determined by instrumental accu- 
racy. Two schemes of the measurements were exploited, 
which correspond to the transmission and back-scattering 
geometries [Fig. 1 (a), the inset]. The measurements 
were complemented by TR and continuous wave (cw) PL 
spectroscopies, as well as cw transmission and reflection 
measurements performed at different temperatures using 
a tungsten lamp. We investigated high-quality c-plane 
ZnO samples with the thickness L of 0.3, 0.4, 1, and 2 
mm, supplied by Tokyo Denpa Co. In the PL spectra, 
the number (2-4) and intensity of the dominant BX lines 
varies among the samples, while their energies were al- 
most identical. 

The typical images of the light pulses propagating at 
different energies are shown in Fig. 1 (a). The leading 
edge of a curved light stripe, gradually narrowing due to 
increasing absorption, is observed up to 3.374 eV with the 
delay about 1.6 ns, i.e. in the close vicinity of the A ex- 
citon. Further shift of the pulse towards higher energies 
results in its full quenching accompanied by the increase 
of BX PL. In the region of relative transparency, the im- 
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FIG. 1: (color online) (a) Combined images of four pulses 
propagating through the 0.3-mm sample (2 K). The bound- 
aries between them are marked by white lines. The insert 
presents the scheme of signal detections, where n • L denote 
the series of replicas, (b) Spectra of cw PL, transmission, and 
reflection (the inset) measured in the same sample. At and 
A l denote the transverse and longitudinal A exciton emission 
peaks. Dashed lines show the transmission spectra simulated 
neglecting the BX lines with hTj\ 1) - 3 /xeV and 2) - 30 /xeV. 



ages contain series of replicas, arising due to the light 
pulse reflection from the crystal boundaries. The mech- 
anism of their transfer is pure ballistic [l7j . They cover 
a distance of n • L, where n is odd (even) number for 
the transmission (back-scattering) geometry. The well- 
defined temporal intervals between them provide a high 
accuracy in determination of the delays. 

The images apparently show that the general curva- 
ture and delay of the pulses follow the optical dispersion 
controlled by the exciton-polariton resonances. The in- 
fluence of the BX lines is local, namely: they provide 
dips cutting the pulses into several parts and induce ex- 
tra light retardation nearby. Because even a very weak 
line can provide noticeable absorption if a sample is thick 
enough, the number of the BX lines resolved by the time- 
of-flight spectroscopy (similar to transmission) is higher 
than in the PL spectra [Fig. 1 (b)]. The spectral cross- 
sections of the pulses (Fig. 2) allow us to reconstruct the 
fine structure of the BX spectrum which contains at least 
18 lines in the 3.356^-3.374 eV range (Table I). In gen- 
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FIG. 2: Selected spectral cross-sections of the light pulses 
passing the 0.3-mm sample, made at different delay times 
from 110 ps (bottom) up to 1200 ps (top). The dips corre- 
sponding to several important lines are labelled. The initial 
light pulse is shown by the dotted curve. 



eral, these lines can be divided into three groups: i) Two 
sets of A and B excitons bound to neutral donors, D°X, 
separated by a gap of ~4.5 meV. ii) Excitons bound to 
ionized donors, D + X, at higher energies, iii) Minor lines 
related to the vibrational-rotational (vr) states of the 
bound excitons. The lines of the first group produce the 
strongest local changes. 

To model the pulse delay and distortion, the initial 
spectrum of the laser pulse E = Eo(uj) in the time mo- 
ment to is approximated as a Gaussian with the cen- 
tral frequency uj$. The amplitude and phase of the pulse 
are given by the Fourier transform at the boundary of a 
medium z = in the linear regime. In accordance with 
the general theory of the electromagnetic field propaga- 
tion, the electric field in the spatial point z = L in the 
moment t can be described as [9[: 



E(L,t) 
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Here, k(u) — (cj/c) y/s{fjj) is the complex wave vector 
of light, e(uj) is the complex dielectric function of the 
medium. 

For the medium with several exciton-polariton reso- 
nances, e(uj) is written as: e(uj) = + Y2j Xj, where £5 
is the background dielectric constant and Xj are the con- 
tributions of these resonances. When the resonances are 
assumed to be homogeneous, the expression taking into 
account the spatial dispersion can be readily used (2o| . 
With inhomogeneous broadening induced, e.g., by struc- 
tural imperfections, each Xj term can be represented by 
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TABLE I: Parameters of bound excitons in bulk ZnO ob- 
tained thought the pulse distortion analysis. The notation 
follows Ref. X denotes unidentified lines. 



Transition 


E (eV) 


fBX 


Transition 


E (eV) 


fBX 


K{D + X A ) 


3.3734 


5 


10~ 7 


X 


3.3635 


5 


10" 7 


h{D+X A ) 


3.3724 


1 


1CT 6 


r Q r {D Q x A ) 


3.3623 


1 


10~ 6 


I 0a (D + X A ) 


3.3718 


1 


1(T 6 


I? r (D°X A ) 


3.3620 


1 


10~ 6 


X 


3.3705 


3 


io- 7 


ir{D°x A ) 


3.3616 


1 


10~ 6 


X 


3.3695 


3 


io- 7 


h(D°X A ) 


3.3614 


2 


10~ 6 


X 


3.3685 


1 


1(T 6 


h(D°X A ) 


3.3608 


5 


10~ 6 


H(D°X B ) 


3.3669 


1 


icn 6 


Ir(D°X A ) 


3.3600 


5 


10~ 6 


lf(D°X B ) 


3.3653 


2 


1(T 6 


h(D°X A ) 


3.3593 


5 


io- 9 


I 7 B (D°X B ) 


3.3647 


2 


1(T 6 


I 9 (D°X A ) 


3.3566 


1 


io- 9 



the convolution of the line with the Gaussian centered on 
the same frequency [l7j |: 

Y _ f fjUQj L_ 

Aj ~ J u j + (3k? + e - iTj - uj 0fa, exp Ui 2 H' 

(2) 

Here, each j resonance is characterized by a frequency 
(jJqj, an oscillator strength /j, and a damping term Tj. 
Aj describes the inhomogeneous width. Spatial disper- 
sion is taken into account by the term f3k 2 = (h 2 /2Mj)k 2 , 
where the effective masses Mj are assumed to be infinite 
for BX and equal 0.9m e for each free exciton resonance. 
Equations (pQ) and ((2]) allow for the spatial dispersion in a 
simplified way, they are valid for uj not too close to ljqj, 
which is sufficient for our purposes. The pulse shapes are 
calculated from (pQ) by means of the Wigner transform. 

In our simulations, all noticeable BX lines were taken 
into account. The HT = 1 jaeY and HA = 75 /ieV charac- 
terize these narrow lines. The Jbx values corresponding 
to the observed local distortion are given in Table I. For 
the sake of illustration, we present the modelling of the 
intricate pattern appearing when the central energy of 
the pulse falls to the gap between the closely situated Iq 
and I7 lines. The light propagates there as an extremely 
narrow ~0.3 meV stripe [Fig. 1 (a)]. Figure 3 shows 
that a twice higher / value taken for the D°X# transi- 
tions (4 • 10 -6 vs 2 • 10 -6 ) provides an unacceptably long 
delay and too strong distortion of the pulse. 

To determine the parameters of the exciton-polariton 
resonances, we consider exclusively the data on the pro- 
cesses taking place inside the crystal, namely: cw trans- 
mission and pulse propagation. We abstain from using 
the data of reflection and PL spectroscopies, because they 
probe mostly a region close to the surface. The mod- 
elling is simplified owing to the fact that each excitonic 
parameter is responsible predominantly for a particular 
characteristics. For instance, the oscillator strength con- 
trols the overall curvature of a pulse via the variation of 
the group velocity v g (uo) = duj/dk and, hence, the dif- 
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FIG. 3: (color online) Simulations of the shape of the 5-ps 
pulse centered at 3.365 eV propagating through the 0.3-mm 
sample: (a) assuming no BX lines; (b,c) taking into account 
the BX lines with (b) / = 2-10 -6 and (c) / = 4-10 -6 . Arrows 
mark the experimental delay of the pulse center. 



ferent delay T = L/v g of the pulse constituents. The 
derived exciton-polariton parameters are collected in Ta- 
ble II (the distant C exciton is neglected). 
The modelling reveals two severe limitations: 

1) The damping term hT must be as small as ~3 fieV. 
Higher values would result in the full opaqueness in the 
transmission spectra in the range where the narrow BX 
lines are clearly resolved [Fig. 1 (b)]. This is indepen- 
dent on the model used to describe e(u)), because the 
absorption depends on the imaginary part of the dielec- 
tric function, which is determined by this term. 

2) The inhomogeneous width HA cannot exceed 0.5 
meV for the A exciton resonances. This restriction arises 
from the observation of the transmitted light at 3.374 eV. 
With stronger broadening, which may occur in crystals 
of worse quality or with a temperature rise, the pulse 
propagating at this energy would be attenuated up to 
full disappearance. The pulse maximum can be at lower 
energy due to enhanced absorption at the leading edge. 

The critical value of the effective damping param- 
eter, yet maintaining the polaritonic modes [2l|, is 
huo^fjhuo/Mc 2 ) 1 / 2 ~ 1.6 meV in ZnO for the A ex- 
citon, being close to the longitudinal-transverse splitting 



TABLE II: Exciton-polariton parameters derived from the 
volume- and surface-probing measurements. 
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FIG. 4: (color online) Light pulses propagating through the 
0.3-mm (a,b) and 1-mm (c,d) samples. The registration is 
done at 2 K in the transmission (a,c) and back-scattering (b,d) 
geometries. The black lines present the delay dependencies 
calculated taking into account the BX transitions. 



ujlt [22]. Our constraints based on the signal attenua- 
tion are more rigorous. In fact, they do not exclude that 
the excitonic resonances inside the crystal are homoge- 
neous with for of a few fieV. With these limitations, we 
have simulated the TR images of the pulses of different 
energies and the series of replicas, using the same exciton- 
polariton parameters for all samples (Fig. 4). 

Our modelling has revealed obvious contradiction be- 
tween the results of volume- and surface-probing exper- 
iments (Table II). The set of the parameters, perfectly 
fitting the cw and TR transmission data, results in too 
sharp peaks in the simulated spectra of reflectivity. The 
successful fitting of the spectrum in Fig. 1 (b), using 
the model of the homogeneous resonances [20], requires 
the empirical damping term hY* A = 0.75 meV. Similar 
result can be obtained assuming the resonances as inho- 
mogeneous with HA a ~ 0.9 meV. However, both variants 
would provide too strong attenuation of the propagating 
light. Further, the splitting between the PL peaks, as- 
cribed to transverse and longitudinal A exciton emission, 
is 1.6 — 2 meV, whereas the rough estimation using the 
bulk parameters gives ujlt ~ 3 meV. There are differ- 
ences in the resonance energies as well. This dependence 
on the measurement technique is suggestive of certain 
broadening and deterioration of the excitonic resonances 
at the surface, which can be induced by modified struc- 



tural properties and local electric fields 

In conclusion, our studies exhibit that a strong reso- 
nant delay (up to 1.6 ns) of light pulses in ZnO is accom- 
panied by their severe attenuation and shape distortion. 
Simulation of the shapes enables us to determine the ex- 
citonic parameters, inherent for bulk ZnO, and to estab- 



lish that the resonances must have the low damping (~ 3 
fjbeV) and limited broadening to allow the light propa- 
gation. We believe that this paper will draw attention 
to the time-of-flight spectroscopy as a new opportunity 
to investigate the fine structure of BX spectrum. The 
discrepancy between the surface and volume parameters 
implies that only the processes taking place inside the 
crystal, like the pulse propagation, are suitable to give 
the true bulk characteristics. 
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